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ABSTRACT 

Magnetic field topology, thermal structure and plasma motions are the three main factors affecting 
the polarization signals used to understand our star. In this theoretical investigation, we focus on the 
effect that gradients in the macroscopic vertical velocity field have on the non-magnetic scattering 
polarization signals, establishing the basis for general cases. We demonstrate that the solar plasma 
velocity gradients have a significant effect on the linear polarization produced by scattering in chromo- 
spheric spectral lines. In particular, we show the impact of velocity gradients on the anisotropy of the 
radiation held and on the ensuing fractional alignment of the Ca II levels, and how they can lead to an 
enhancement of the zero-field linear polarization signals. This investigation remarks the importance 
of knowing the dynamical state of the solar atmosphere in order to correctly interpret spectropolari- 
metric measurements, which is important, among other things, for establishing a suitable zero field 
reference case to infer magnetic fields via the Hanle effect. 
Subject headings: Polarization - scattering - radiative transfer - Sun: chromosphere, quiet sun 



1. INTRODUCTION 

Over the last few years it has become increasingly 
clear that the determination of the magnetic field in the 
"quiet" solar chromosphere requires measuring and in- 
terpreting the linear polarization profiles produced by 
scattering in strong spectral lines, such as H Q and the 
8542 A line of the infrared (IR) triplet of Ca II (e.g., sec 
reviews by Trujillo Bueno 2010 and Uitenbrock 2011). In 
these chromosphcric lines, the maximum fractional lin- 
ear polarization signal occurs at the center of the spectral 
line under consideration, where the Hanle effect (i.e., the 
magnctic-ficld-induccd modification of the scattering line 
polarization) operates (Stenflo 1998). Since the opacity 
at the center of such chromospheric lines is very signif- 
icant, it is natural to find that the response function 
of the emergent scattering polarization to magnetic field 
perturbations peaks in the upper chromosphere (Stcpan 
Sz Trujillo Bueno). This contrasts with the circular po- 
larization signal caused by the Zeeman effect whose re- 
sponse function peaks at significantly lower atmospheric 
heights (Uitcnbroek & Socas Navarro 2004) . Of particu- 
lar importance for developing the Hanle effect as a diag- 
nostic tool of chromospheric magnetism is to understand 
and calculate reliably the linear polarization profiles cor- 
responding to the zero-field reference case. 

The physical origin of the scattering line polarization 
is atomic level polarization (that is, population imbal- 
ances and/or coherence between the magnetic sublevels 
of a degenerate level with total angular momentum J). 
Atomic polarization, in turn, is induced by anisotropic 
radiation pumping, which can be particularly efficient in 
the low-density regions of stellar atmospheres where the 
depolarizing role of elastic collisions tends to be negli- 
gible. The larger the anisotropy of the incident field, 



ccarlin@iac.es 

1 Institute- de Astrofi'sica de Canarias, 38205, La Laguna, 
Tenerife, Spain 

2 Departamento de Astrofisica, Universidad de La Laguna, 
Tenerife, Spain 

3 Consejo Superior de Investigaciones Cientificas, Spain 



the larger the induced atomic level polarization, and the 
larger the amplitude of the emergent linear polarization. 
The degree of anisotropy of the spectral line radiation 
within the solar atmosphere depends on the center-to- 
limb variation (CLV) of the incident intensity. In a static 
model atmosphere the CLV of the incident intensity is 
established by the gradient of the source function of the 
spectral line under consideration (iTruiillo Bueno! [20011 : 
ILandi Pegl'Innocenti fc Lan dolfi 2004|). However, stellar 
chromospheres are highly dynamic systems, with shocks 
and wave motions (e.g., Carlsson & Stein 1997). The en- 
suing macroscopic velocity gradients and Dopplcr shifts 
might have a significant impact on the radiation field 
anisotropy and, consequently, on the emergent polariza- 
tion profiles. Therefore, it is important to investigate the 
extent to which macroscopic velocity gradients may mod- 
ify the amplitude and shape of the emergent linear po- 
larization profiles produced by optically pumped atoms 
in the solar atmosphere. The main aim of this first pa- 
per is to explain why atmospheric velocity gradients may 
modify the anisotropy of the spectral line radiation and, 
therefore, the emergent scattering line polarization. We 
aim also at evaluating, with the help of ad-hoc velocity 
fields introduced in a semi-empirical solar model atmo- 
sphere, their possible impact on the scattering polariza- 
tion of the IR triplet of Ca II. A recent investigation by 
Manso Sainz & Trujillo Bueno (2010), based on radiative 
transfer calculations in static model atmospheres, shows 
why the differential Hanle effect in these lines is of great 
potential interest for the exploration of chromospheric 
magnetism. 

2. FORMULATION OF THE PROBLEM AND RELEVANT 
EQUATIONS 

2.1. The atomic model and the statistical equilibrium 
equations (SEE) 

We assume an atomic model consisting of the five low- 
est energetic, fine structure levels of Ca II (see Figure 
1). The excitation state of the atomic system is given 
by the populations of its 18 magnetic sublevels and the 
coherences among them. We neglect coherences between 
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TABLE 1 
Short list of Ca ii atomic data parameters. From left to 
right: the central wavelength, the upper and the lower 
level of each transition, the radiative rates from nist 

atomic spectra database 

(http://www.nist.gov/physlab/data/asd.cfm) AND THE ATOMIC 

polarizability coefficients introduced by landi 

Degl'Innocenti (1984). 



different energy levels of the same term (multilevel ap- 
proximation). Moreover, since the problem we consider 
here (plane-parallel, non-magnetic atmosphere with ver- 
tical velocity fields) is axially symmetric around the ver- 
tical direction, no coherences between different magnetic 
sublcvels exist when the symmetry axis is taken for quan- 
tizing the angular momentum. We use the multipolar 
components of each J-lcvel, 
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(1) 
where K = 0,...,2J, Nm is the population of the sub- 
level with magnetic quantum number M , and the symbol 
between brackets is the Wigner 3j-symbol (e.g., Brink & 
Satchler 1968). Due to the symmetry of the scattering 
process (no magnetic field, no polarized incident radi- 
ation in the atmosphere's boundaries), in a given level 
N + m = N-m, and the excitation state of the system 
is described by just 9 independent sublevel populations. 
Consequently, odd-K elements (orientation components) 
in Equation 1 vanish for all five levels, and the only inde- 
pendent variables of the problem in the spherical compo- 
nents formalism are the total populations of the five lev- 
els (y/2J + lpo)i th° alignment components (p 2 ,) of levels 
2,3, and 5; and p$ of level 3, whose role is negligible for 
our problem. 

The statistical equilibrium equations accounting for 
the radiative and collisional excitations and deexcita- 
tions in the 5-level system of Figure 1 are given explic- 
itly in Manso Sainz & Trujillo Bueno (2010 ;hereafter 
MSTB2010). We have particularized them to the no- 
coherence case (only p$ elements) in Appendix [B] The 
statistical equilibrium equations for the p® components 
contain terms that are equal to those appearing in the 
statistical equilibrium equations for the populations in a 
standard (no polarization) NLTE problem (e.g., Mihalas 
1978), plus higher order terms ~ Jq/9 2 (see Eqs, IBlllB5|) . 
The statistical equilibrium equations for the alignment 
(Pq components) are formally similar to the ones for the 
populations with additional terms ~ J0P0 accounting for 
the generation of alignment from the anisotropy of the ra- 
diation field, and (negligible) higher order terms ~ JqPq 
and J%pl (see Equations (|B6]| . (|B7|) and (jBSj)). These 
equations are expressed in the atom reference frame (co- 
moving system). 

Since the radiation field is axially symmetric, just two 
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Fig. 1. — Atomic model with energy levels for Call. The labels 
indicate the spectroscopic terms. Numbers inside blue filled circles 
identificate each level. Solid lines connecting levels show the al- 
lowed radiative transitions and the numbers at the middle of each 
segment are their wavelength in A. The wider the width of each 
connecting line the larger the spontaneous emission rate A u ; of the 
tran sition. Atomic data for each spectral line are shown in Table 

i2~n 

radiation field tensor elements ( Jq and Jq) are necessary 
to describe the symmetry properties of the spectral line 
radiation. Let I(v,p) and Q(v,p) be the Stokes param- 
eters expressed in the observer's frame at a given height 
z, where v is the frequency, p = cos 9 and 9 is the an- 
gle that the ray forms with the vertical direction. Then, 
the corresponding values seen by a comoving frame with 
vertical velocity v z with respect to the observer's frame 
arc I'{v',p) = I{v,p) and Q'{v',p) = Q(i/,p), where 
1/ = v{l—v z p/c) and v = v'(l + v z p,/c) (to first order in 
v z /c). Therefore, the mean intensity at the considered 
height, can be expressed from one or another reference 
frame as: 
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where 4>tu{v) is the absorption profile (e.g., for 
a Gaussian profile, we would have 4>t u {y) = 
7t -1 / 2 Ai/d -1 exp(— (y — ^ ) 2 /Ai/£) 2 ), with v Q the cen- 
tral line frequency and Az/q the Dopplcr width) and 
4>' lu (v,p) = 4>( u (v{\ - v z p/c)), with v z > for upflowing 
material. Analogously, the anisotropy in the observer's 
frame: 
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The important quantity that controls the ability of 
an anisotropic radiation field to generate atomic polar- 
ization is the line anisotropy factor for each transition, 
which can be calculated as: 
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Its range goes from iyu ne = —0.5 (for a radiation field 
coming entirely from the horizontal plane) to wii nc = 1 



(for a collimated vertical beam). 

2.2. The radiative transfer equations (RTE) 

Due to symmetry, in a non-magnetized plane-parallel 
medium with a vertical velocity field, light can only be 
linearly polarized parallel or perpendicularly to the stel- 
lar limb. Therefore, chosing the reference direction for 
positive Q parallel to the limb, the only non-vanishing 
Stokes parameters are I and Q, and they satisfy the fol- 
lowing radiative transfer equations: 
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where s is the distance along the ray. The absorption 
and emission coefficients are (MSTB2010): 
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where r/i c ° nt and e 7 cont are the continuum absortion and 
emission coefficients for intensity, respectively Likewise, 
rjj e and t/q 10 are the line absorption coefficients for 

Stokes / and Q, respectively, while e' 7 lnc and e^ 10 are 

the line emission coefficients for Stokes / and Q, respec- 

(2) (2) 

tively. The coefficients Wj j and Wj^j depend only on 

the transition and are detailed in Table (|2.1j) .The sub- 
scripts u and £ refer to the upper and lower level of the 
transition considered, respectively, and B v is the Planck 
function at the central frequency v$ of the transition. 
Note also that: 
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where J\f is the total number of atoms per unit volume. 
With the total absortion coefficient for the intensity, 
the line of sight (los) optical depth for each frequency is 
calculated by the following integral along the ray: 
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2.3. Numerical method 

The solution to the non-LTE problem of the sec- 
ond kind considered here (the self-consistent solution 
of the statistical equilibrium equations for the den- 
sity matrix elements together with the radiative trans- 
fer equations for the Stokes parameters) is carried out 
by generalizing the computer prog ram developed by 
iManso Sainz fc Truiillo Buenol (|2003| ) , to allow for radial 
macroscopic velocity fields. For integra ting the RTE, a 
parab olic short-characteristics scheme (|Kunasz fc Auerl 
119881 ) is used. At each iterative step, the radiative trans- 
fer equation is solved, and Jq and Jj are computed 
and used to solve the SEE following the accelaratcd 
Lambda iteration method outlined in the Appendix of 
MSTB2010. Once the solution for the multipolar com- 
ponents of the density matrix is consistently reached, the 
emergent Stokes parameters are calculated for the de- 
sired line of sight, which in all the figures of this paper is 
fi = 0.1. This final step is done increasing the frequency 
grid resolution to a large value in order to correctly sam- 
ple the small features and peaks of the emergent profiles. 

Some technical considerations have to be kept in mind 
for the treatment of velocity fields. Due to the presence 
of Doppler shifts, the wavelength axis used to compute 
Jq and Jq must include the required extension and res- 
olution, because the spectral line radiation may be now 
shifted and asymmetric. In our strategy for the wave- 
length grid, the resolution is larger in the core than in the 
wings, keeping the same wavelength grid for all heights. 

The cutoff wavelength for the core (where resolution 
is appreciably higher) is dictated by the maximum ex- 
pected Doppler shift. Thus, the core bandwidth is es- 
timated allowing for a range of 2V, nax around the zero 
velocity central wavelength of the lines, with Vmax the 
maximum velocity found in the atmosphere (in Doppler 
units). Apart from that, a minimum typical resolution 
for the core is set to 2 points per Doppler width (Ai/p). 
Then, the height with the smallest Doppler width de- 
termines the core resolution, and the height with the 
maximum macro-velocity states its bandwidth. 

Furthermore, as frequencies and angles are inextrica- 
bly entangled (through terms v — v z \xv /c appearing in 
the absortion/emission profiles due to the Doppler ef- 
fect, like in Eq. ((21)), the maximum angular increment 
(A/^ max ) is restricted by the maximum frequency incre- 
ment (Ax max w 1/2, in Doppler units). Thus, it must oc- 
cur that A^t max • Vmax < 1/2. In the worst case, the max- 
imum allowed angular increment will be smaller (more 
angular resolution needed) when the maximum vertical 
velocity increases. Besides this consideration, the maxi- 
mum angular increment could be even more demanding 
because of the high sensitivity of the polarization profiles 
to the angular discretization. 

Finally, the depth grid must be fine enough, in such a 
way that the maximum difference in velocity between 
consecutive points is not too large, the typical differ- 
ence being equal to half the Doppler width (|V(zj) — 
V"(«i-i)| < 1/2). If the difference is larger, the ab- 
sortion/emission profiles would change abruptly with 
height, producing imprccisions in the optical depth in- 
crements (Mihalas, 1978). 
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3. EFFECT OF A VELOCITY GRADIENT ON THE 
RADIATION FIELD ANISOTROPY AND THE MEAN 
INTENSITY 

As we shall see, the presence of a vertical velocity gra- 
dient in an atmosphere enhances the anisotropy of the 
radiation field and, hence, of the scattering line polariza- 
tion patterns. The fundamental process underlying this 
mechanism can be simply understood with the following 
basic examples. 

3.1. Anisotropy seen by a moving scatterer. 

Consider an absorption spectral line with a Gaussian 
profile emerging from a static atmosphere with a linear 
limb darkening law, 



I(v,H) = I (0) (1 - u + un)[l -oexp(-(- 
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where I^ ' is the continuum intensity at disk center, u 
is the limb darkening coefficient, a < 1 measures the 
intensity depression of the line, and w its width. In 
this approximation we assume that all the parameters 
are constant. Now imagine that, at the top of the at- 
mosphere, there is a thin cloud scattering the incident 
light given by Eq. ([10[) and moving radially at velocity 
v z with respect to the bottom layers of the atmosphere, 
supposed static. We will assume that the absorption pro- 
file is Gaussian (dominated by Doppler broadening) , with 
width Aj/£>. When Ai^rj <C w, the incident spectral line 
radiation is much broader than the absorption profile (in 
fact, for a = A.i/£>/w = the absorption profile is for- 
mally a Dirac-(S function). Then, from Equations (O-jSJ) 
we can derive explicit expressions for the mean inten- 
sity and anisotropy of the radiation field as a function of 
the scatterer velocity (see Appendix A): Jq = Xo(a; £)/2, 
Jo = 22(c*;£)/4\/2, where the Xo,2 functions are defined 
by Equations (|A"4 ]) -(|A"5 j) . The behaviour of Jq and Jq 
with the adimensional velocity £ = v z vqJ{cw) (Figure 
[2]), is most clearly illustrated in their asymptotic limits 
at low velocities: 
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Equation (|11|) shows that, for an absorption line (a > 0), 
Jq is always increasing with the velocity since the co- 
efficient of £ 2 is positive, and this, regardless of the 
sign of v z (i.e., regardless of whether the scatterers 
move upwards or downwards); if the line is in emission 
(a < 0), Jq monotonically decreases. These are the 
Doppler brightening and Doppler dimm ing effects (e.g., 
ILandi Degl'Innocenti fc Landolfil I2004D . An analogous 
analysis applies to Jq (Eq. [T2"|) . Note that, in the ab- 
sence of limb darkening (u = 0), the anisotropy vanishes 
in a static atmosphere, while the mere presence of a rel- 
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Fig. 2. — J§//(°) (upper panel) and \/2J§/J§ (lower panel) as 
a function of the adimensio nal velocity £ calculated using an in- 
cident line profile as in Eq. I|10j l with u = 0.3 and a = 0.4 or 0.8 
(see labels). Dotted lines have been computed for the case of an 
infinitely sharp absorption profile (a = 0). Solid lines refer to the 
case a = 0.9 (non saturated line). The case with no limb darken- 
ing (u = 0) and a = 0.4 has been plotted for comparison (dashed 
lines). 

ative velocity between the scatterers and the underly- 
ing static atmosphere induces anisotropy in the radiation 
field — hence, a polarization signal. A real atmosphere 
could then be understood as a superposition of scatterers 
that modify the anisotropy depending on the local ve- 
locity gradient and the illumination received from lower 
shells. An interesting discussion on the effect of velocities 
with directions other than radial can be found in Landi 
Degl'Innocenti & Landolfi (2004; Section 12.4). 

Clearly, all the above discussion depends on the 
Doppler shift induced by the velocity v z normalized to 
the width of the spectral line, i.e., on £. A large veloc- 
ity gradient on a broad line has the same effect as that 
of a smaller velocity gradient on a narrow line. This is 
important to be kept in mind since the response of dif- 
ferent spectral lines to the same velocity gradient will 
be different, what can help us to decipher the velocity 
stratification. Even different spectral lines belonging to 



the same atomic species may have very different widths, 
as for example, the Ca n IR triplet and the UV doublet 
studied in the next section. 

3.2. Calculations in a Milne-Eddington model. 

The discussion above explains the basic mechanism 
by means of which a velocity gradient enhances the 
anisotropy of the radiation field. Now we can get fur- 
ther insight on the structure of the radiation field within 
an atmosphere with velocity gradients from just the for- 
mal solution of the RT equation for the intensity (e.g., 
Mihalas 1978). As before, we neglect effects due to po- 
larization and Jq and Jq are calculated from Stokes / 
alone. We consider a semi-infinite, plane-parallel atmo- 
sphere with a source function S = Sq(1 + /3ti), where 
77 is the integrated line optical depth in the static limit 
(hence, the element of optical depth dr v — (r + <f>[v(l — 
v z (Ti)n/c)])dTi, where r — k c /ki is the ratio of contin- 
uum to line opacity). We begin by considering a ver- 
tical velocity field v z (ti) = v /[l + (ti/t )] (positive 
outward the star), shown in Fig. [3l Equivalently, we 
may express the velocity in adimensional terms by using 
£ = (vo/cjVz/Aud (the width Avd of the Gaussian ab- 
sorption profile is assumed to be constant with depth). 
The parameter To fixes the position of the maximum ve- 
locity gradient. Note that the wavelength dependence 
of the Doppler effect (AX Z = \qv z /c) is cancelled in the 
adimensional problem, where velocities are measured in 
Doppler units. It is easy to calculate numerically I g (y, /i) 
at every point in the atmosphere and thus, the mean in- 
tensity and anisotropy of the radiation field (Figure 2]). 

The rise in Jq in higher layers with respect to the static 
case corresponds to the Doppler brightening discussed 
above. Thanks to the Doppler shifts, the atoms see more 
and more of the brighter continuum below, which en- 
hances Jq. When the maximum velocity gradient takes 

place at optically thick enough layers (to > 1), Jq is also 
larger than for the static case, but it decreases monoton- 
ically with height in the atmosphere (to = 10 2 , dotted 
lines in upper panel of Fig. [4|). Note that the impor- 
tant quantity that modulates the increase in JS is not 
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Fig. 3. — Adimensional velocity fields considered in this sec- 
tion. They are parameterized by the limiting value at small optical 
depths (labels) and by the location in optical depth of the largest 
velocity gradient region (vertical lines marking to). Solid lines: 
to = 1. Dotted lines: to = 100. 
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Fig. 4. — Jq (upper panel) and V^Jq/Jq (lower panel) as a 
function of the integrated static line opacity (tj ) in an expanding 
atmosphere with S = S^°'(l + firi) and different velocity stratifica- 
tions £ = $o/(l + tj/to) (see Fig. [5). The parameters used in this 
plot are /3 = 3/2 and k c /ki = 10 -4 . The labels indicate the value 
of £o- Solid lines correspond to to = 1 while dotted lines refer to 
to = 100. 

the maximum velocity but the velocity gradient (differ- 
ence in velocity between optically thick and optically thin 
parts of the atmosphere). The larger the gradient, the 
more pronounced the radiative decoupling is between dif- 
ferent heights. An extreme example of such radiative de- 
coupling could be found in supernovae explosions, where 
the vertical velocity gradients are huge. 

In our case (vertical motions), the Doppler brighten- 
ing implies an enhancement of the contribution of verti- 
cal radiation to Eq. <j3j) with respect to the horizontal 
radiation, with the latter remaining almost equal to the 
static case (no horizontal motions, no horizontal Doppler 
brightening) . This velocity- induced limb darkening is the 
origin of the anisotropy enhancement. 

However, note that the maximum anisotropy does not 
rise indefinitely when increasing the maximum velocity. 
If the velocity gradient in units of the Doppler width is 
larger than ~ 3 (see curves for £ max = 5 in Fig. [4j, the 
anisotropy at the surface saturates and decreases (even 
below the curves corresponding to shorter velocity gra- 
dients). It is accompanied by a bump around 77 = 1 
when the maximum velocity gradient is taking place at 
low density layers (to < 1). This behavior can be under- 
stood using Eq. ^. When ( max < 3, an increment in 
Cmax entails a rise in Jq , Jq and Jq/ Jq (w\i nc ) in the up- 
per atmosphere, what means that the velocity gradients 
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Fig. 5. — Same as Fig. [4] for an isothermal moving atmosphere 
with a two-level atom model using the velocity fields of Fig. [3] We 
assume a very strong line (r c = 0) and e = 10 — 4 . Upper panel: 
the Doppler brightening effect increases the surface source function 
value. In a two-level model atom in a static atmosphere (label 0) 
this value follows the well known expression S(0) = \fkB. The 
vertical axis is in units of the Planck function. Lower panel: we 
show the amplification of the radiation field anisotropy when the 
velocity gradient increases. The vertical line marks the location of 
the largest velocity gradient (r; = To = 0.1). The curve labeled 
with "5" shows the saturation of the anisotropy and the bump 
produced by the strong velocity gradient taking place at optically 
thin regions. 

enhance the imbalance between vertical and horizontal 
radiation. However, if £ max is above that threshold, Jq 
and Jq rise, but the ratio Jq/ Jq saturates and dimin- 
ishes. The reason is that a large velocity gradient makes 
the absorption profiles associated to almost horizontal 
outgoing rays (0 < fj, < l/\/3) shifted so much that they 
also capture the background continuum radiation. Then 
contributions are negative to the angular integral of Jq 
but positive for Jq. 

Separating the contributions of rays with angles in the 
range l/\/3 < |/i| < 1 (that we refer to with the label +) 
and angles in the range < |/i| < l/y3 (that we refer 
to with the label — ), the line anisotropy can be written 

aS Wlinc = ' 
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appreciably when £ max > 3. Therefore, although iyii ne 
increases for all velocity gradients, its enhancement is 
smaller for large velocity gradients than for smaller ones. 
This effect occurs as well when motions take place deeper 
(to > 1) but it is less important and the anisotropy bump 
and saturation arc reduced. 



For the considered velocity fields (with a negligible 
gradient in the upper atmosphere), Jq and Jq reach an 
asymptotic value in optically thin regions. We have veri- 
fied that this effect does not occur if the velocity gradient 
is not zero in those layers. In any case, the presence of 
a large anisotropy in optically thin heights barely affects 
the emergent linear polarization profiles. 

3.3. Two-level model atom in dynamic atmospheres. 

Before going to a more realistic case, a final illustrative 
example is considered. In this case, we assume the same 
parameterization of the velocities than in the previous ex- 
ample, but now we solve the complete iterative RT prob- 
lem with a two-level atom model and a specific temper- 
ature stratification. Consequently, the source function 
and the anisotropy are consistently obtained in a mov- 
ing atmosphere. The intens ity source function is Sj 



Sf nc + (1 - r Ufl )B (e.g., IRvbicki fc Hummed U992T) . 



with r„ M = <p' eu (iy,fj)/(r c + <f>' iu (v, fi)) and the expression 
for the line source function remains formally equal to that 
of the static case, being 5} mc = (1 — e) Jq +eB, where B is 
the imposed Planck function, e is the inelastic collisional 
parameter and Jq is calculated with Eq. (J2J|. 

The qualitative behaviour explained in the previous 
subsection is maintained in these two-level atom calcula- 
tions. For small e values (large NLTE effects), the source 
function Sj ~ Jq shows Doppler brightening effects and 
its surface value depends on the maximum velocity gra- 
dient and on the maximum background continuum set 
by the photospheric conditions (upper panel in Fig. [5|. 
The anisotropy rises proportionally to the velocity gra- 
dient until a saturation occurs (lower panel in Fig. [5]). 
A similar behavior is found when the maximum velocity 
gradient occurs higher in the atmosphere (sec Fig. 1101 in 
Appendix [C| . 

In a static atmosphere, the radiation field anisotropy 
is dominated by the presence of gradients in the in- 
tensity source function (Trujillo Bueno 2001; Landi 
Degl'Innocenti & Landolfi 2004), which can be modi- 
fied via the Planck function (equivalently, the temper- 
ature). In the dynamical case that we are dealing with, 
the slope of the source function is also modified due to the 
existence of velocity gradients thanks to the frequency- 
decoupling caused by relative motions between absortion 
profiles (Doppler brightening). In general, both mech- 
anisms act together (velocity-induced and temperature- 
induced modification of the source function gradient) and 
the ensuing anisotropy and the emergent linear polariza- 
tion profiles are modified accordingly. 

It is important to note that the adimensional velocity £ 
depends both on the velocity and also on the line Doppler 
width because £(t) = 5v / Ai/p = v z / \f2ksT /m, with ks 
the Boltzmann constant, T the temperature and m the 
mass of the atom. In the photosphere, where velocities 
are much lower than in the chomosphere, £ is expected 
to be negligible. In the chromosphere, plasma motions 
are important and the temperature is still comparable to 
that of the photosphere, inducing £ to be controlled by 
the velocity field. However, for layers in the transition 
region and above, the high temperatures reduce the value 
of £. In any case, at these heights, the density is so low 
that, although £ (and consequently the anisotropy) could 
have a highly variable behaviour, the emergent polariza- 
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Fig. 6. — Amplification of the line anisotropy (wi; nc = V2Jq/Jq) due to vertical velocity gradients. Upper left panel: linear velocity fields 
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of t 1 ™ = 1 for the transitions 8498A (blue), 8542A (green), 8662A (red) and the K line (black). Remaining panels: corresponding line 
anisotropy factors plotted against t'° s for each line. 



tion profiles of chromosphoric lines will not be sensitive 
to them. 

4. RESULTS FOR THE CA II IR TRIPLET . 

Now, we study the effect of the velocity field on a 
multilevel atomic system in a realistic atmospheric 
model, within the more general framework described 
in Sec. [2j We consider the formation of the scattering 
polarization pattern of the Ca II infrared triplet in 
a semi-empirical model atmosphere (FAL-C model of 
Fontenla, Avrett, and Loeser 1993) in the presence 
of vertical velocity fields (v = v z (z)k, with k the 
unit vector along the vertical pointing upwards). We 
will assume a constant microturbulcnt velocity field 
of 3.5 km s _1 , a representative value for the region 
of formation which gives a realistic broadening of the 
triplet profiles. 



4.1. Behaviour of the anisotropy in the Ca II IR triplet 

For simplicity, we set linear velocity fields (constant 
velocity gradient along z) between z = —100 and 2150 
km (see Fig. [5]). Consequently, the adimensional velocity 
field £ z has a non-monotonic behavior due to its depen- 
dence on the temperature (upper right panel in Fig. [5]). 
In the chromosphere, where the Ca II triplet lines form, 
£ z is dominated by the macroscopic motions. Here, the 
velocity gradients produce variations in the anisotropy of 
the triplet lines that agree with the behaviour outlined 
in the previous sections. Namely, an amplification and 
a subsequent saturation of the anisotropy factor due to 
the significant velocity gradient at those heights (see the 
lower panels and middle right panel in Fig. [6]). Above the 
chromosphere, on the contrary, the temperature domi- 
nates (£z stabilizes and diminishes) and the anisotropy 
slightly decreases with height. 

If the (adimensional) velocity gradient is negligible 



where the line forms (around r, 



los 



1), the anisotropy 



remains unaffected. Otherwise, if a spectral line forms 
at very hot layers, where the absortion profiles are wider 
and their sensitivity to the velocity gradients is lower, the 
Doppler brightening will not be so efficient amplifying the 
anisotropy. This is the case of the anisotropy of the Ca 
II K line (middle left panel in Fig. [6J. Compare how the 
slope of £ z is smaller where the Ca n K line forms (black 
line on Fig. [6]) than where the triplet lines do. Conse- 
quently, the enhancement of the line anisotropy through 
the presence of velocity gradients in this line is reduced. 
All our calculations demonstrate that the anisotropy 
in the Ca II IR triplet can be amplified through chromo- 
sphcric vertical velocity gradients. This results suggest 
that the same occurs with the ensuing linear polarization 
profiles. 

4.2. The impact on the polarization of the emergent 
radiation. 

For investigating the effect of vertical velocity fields on 
the emergent fractional polarization profiles we perform 
the following numerical experiments. First, we impose 
velocity gradients with the same absolute value but op- 
posite signs (top left panel of Fig. [7]). The resulting 
emergent Q/I profiles (remaining left panels of Fig. [7J 
are magnified by a significant factor (> 2 for all the tran- 
sitions) with respect to the static case (black dotted line). 



The linear polarization profiles have the same amplitude, 
independently of the sign of the velocity gradient. An- 
other remarkable feature is the asymmetry of the profiles, 
having a higher blue wing in those cases in which the ve- 
locity gradient is positive and a higher red wing when 
the velocity gradient is negative, independently of the 
velocity sign. Note also that the Q/I profile is shifted in 
frequency due to the relative velocity between the plasma 
and the observer. 

As a second experiment, we consider different velocity 
fields with increasing gradients (right upper panel in Fig. 
UJ. In the ensuing Q/I profiles we see that the larger the 
velocity gradient, the larger the frequency shift of the 
emergent profiles and the larger the amplitude. In all 
transitions, one of the lateral lobes of the signal remains 
almost constant. Thus, what really changes is the central 
part of the profiles, being a "valley" in the A8498 line and 
a "peak" in the other two transitions. To quantify these 
variations, we define {Q/I) pp (peak-to-peak amplitude of 
Q/I) as the difference between the lowest and the highest 
value of the emergent Q/I signal, which is also a measure 
of its contrast. Note that, as expected from the first 
experiment, (Q/I) pp depends only on the absolute value 
of the gradient. Figure [5] summarizes these results. 

The sensitivity of the linear polarization to the ve- 
locity gradient can be measured approximately as com- 
mented in Sec. 13.11 using a parameter a = Avd/w 
that accounts for the difference in width of the absorp- 
tion profile with respect to the emergent intensity pro- 
file. If a ~ 1, small adimensional velocities will produce 
large changes in shape; if a <C 1, much larger £ z values 
are needed for the same effect. In the case of the IR 
triplet lines, a ~ 0.355 in the formation region of A8498 
and around 0.29 and 0.285 in the formation region of 
A8542 and A8662 (having wider profiles), respectively. 
Then, the former is more sensitive to velocity variations 
in its formation region (Fig. [5]). Finally, the K line has 
a (tI° s = 1) ~ 0.015, a low value due to its wider spectral 
wings. 

The enhancement of the polarization signals are a con- 
sequence of the increase in the anisotropy. Therefore, 
since this increase is produced by the presence of velocity 
fields, the polarization signals of the Ca n IR triplet are 
sensitive also to the dynamic state of the chromosphere. 

4.3. Variations on the atomic alignment due to velocity 

gradients. 

In order to get physical insight on the formation of the 
emergent polarization profiles, we use an analytical ap- 
proximation. Following Trujillo Bucno (2003), the emer- 
gent fractional linear polarization for a strong line at the 
central wavelength can be approximated with: 



2\/2 



(l-^ 2 ) 



■w 



(2) 



■Oq(J u ) 



w 



(2) 
JeJv. 



°a(Jt) 



,(2) 



(13) 
The symbols w^jj, are numerical coefficients already in- 
troduced in Sec. [H The quantities cr 2 (J„) and (Jq{Ji ) are 

2 - pi/pi) eval- 



the fractional alignment coefficients (<7g 



uated at r, 



los 



1 for the upper and lower level of the 



transition, respectively This is the generalization of the 
Eddington-Barbicr (EB) aproximation to the scattering 
polarization and establishes that changes in linear polar- 
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Fig. 7. — Left panels: calculation at fi = 0.1 of the emergent Q/I polarization signals of the Ca II IR triplet when four different choices 
for the vertical velocity gradients are imposed. Positive velocities imply upflowing plasma. Gradient "a"/"b" simulates an atmosphere 
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Each curve is computed on the converged solution of the multilevel NLTE problem described in Sec. [21 Right panels: same calculations 
than in the left-hand panels, but with different velocity gradient values varying from to 16.3ms _1 km -1 in steps of 2.3ms" 1 km -1 (sec 
top right panel). These results show that the polarization signals are increased and shifted with respect to the static case depending only 
on the absolute value of the vertical velocity field gradient and independently of the sign of the velocity field. 



ization (for a static case) are induced by changes in the 
atomic aligment of the energy levels. 

Our calculations show that vertical velocity fields with 
moderate gradients (< 10ms _1 km~ 1 in a linear velocity 
field, as the ones shown in the figures) do indeed produce 
variations in the fractional alignment, which are small 
for |f7o(J u )| and significant for |ctq(J^)| (see Fig. [9]). The 



lower level alignment is the main driver of the changes 
produced in the Q/I signals. This is strictly true for the 
A8662 line, whose upper level with J = 1/2 cannot be 
aligned (zero field dichroism polarization). In the other 
transitions, a certain influence of the upper level align- 
ment becomes important only for large gradients. The 
reason is that the K transition is so strong in comparison 
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Fig. 8. — {Q/I)pp normalized to the static case solution versus 
the velocity gradient (bottom axis) and the maximum absolute 
velocity (top axis) for linear velocity fields appearing in Fig. The 
results are invariant under velocity sign changes. The transition 
A8498 is more sensitive to velocity variations due to its thinner 
absortion profile. 



with the IR triplet that it is dictating the common upper 
level 5 alignment (Fig. [1]). In fact, ctq{J^) is driven by 
the K line anisotropy which, at cromosphcric heights, is 
almost unaffected for the considered velocity gradients, 
as we discussed in Sec. 14.11 (Fig. [6|). Thus, the strong 
H and K lines feed population to the upper levels and 
the K line controls the alignment of the 2 P^/2 level (see 
Fig. 1), while the polarization signals of the IR triplet 
change with velocity fields affecting Oq(J() (through the 
anisotropy enhancement). 

To illustrate the well-known link between the align- 
ment and the anisotropy we can follow the next reason- 
ing. For the Ca n model atom that we deal with in this 
work, it is posible to derive a simple analytic expression 
that relates the anisotropy and the alig nmen ts for the 
A8542 transition. Making use of Eqs. (|lp and (|B5JI 
and neglecting second order terms and collisions, we find 
that: 

2<^(J 5 ) - V7<J 2 (J 3 ) * W Une (3 -> 5). (14) 

As before, we can roughly assume that ctq(Js) ~ constant 
at the chromosphere because it is controled by the K line. 
Then, Eq. (fT4|) suggests that, if the radiation anisotropy 
increases at that heights, an amplification of |<Tq(J 3 )| oc- 
curs (note that cr 2 {Jz) is negative for these lines). A 
more aligned atomic population produce a more intense 
scattering polarization signal. 

5. CONCLUSIONS 

When vertical velocity gradients exist, the polarization 
profiles are always shifted in wavelength, asymmetrized 
and enhanced in amplitude with respect to the con- 
stant velocity case. The reason is that increments in 
the absolute value of the velocity gradient increase the 
source function (Dopplcr brightening) and enhance the 
anisotropy of the radiation field (Sees. [3]and[3|), that 
in turn modify the fractional alignment (Sec. I4.3[) and 
amplify the scattering polarization profiles (Sec. 14 21) . 



For this very reason, all calculations assuming static 
models in the formation region might underestimate the 
scattering polarization amplitudes and not capture the 
right shape of the profiles. In particular, it must be 
taken into account that the Ca n IR triplet lines form un- 
der non-LTE conditions in chromosphcric regions, where 
velocity gradient may be significantly large due to the 
upward propagation of waves in a vertically stratified 
atmosphere (e.g., Carlsson & Stein 1997). Probably, 
in photospheric and transition region lines the effect of 
velocities on polarization can be safely neglected (they 
will be predominantly amplified by temperature gradi- 
ents as discussed in Sec. 13. 3p . but not necessarily in 
the chromosphere. In our study we see that the A8498 
line is more sensitive to macroscopic motions in the low- 
chromosphere, while the A8542 and A8662 lines are espe- 
cially amplified when strong velocity gradients are found 
at heights around 1.5 Mm and higher in our model. 

At the light of these results, it is obvious that the effect 
of the velocity might be of relevance for measuring chro- 
mosphcric magnetic fields. In particular, the described 
mechanism might turn out to be important for the cor- 
rect interpretation of polarization signals in the Sun with 
the Hanle effect. Given that weak chromosphcric mag- 
netic fields are inferred with the Hanle effect using the 
difference between the observed linear polarization signal 
and the signal that would be produced in the absence of 
a magnetic field, it is crucial to correctly compute the 
reference no-magnetic signal. In the Ca n IR triplet, it 
could be possible to break the degeneracy of the com- 
bined effects by taking into account the different sensi- 
tivities of the three lines of the triplet to the magnetic 
field and the thermodynamics. As stated by MSTB2010, 
the A8498 line is very sensitive to the thermal structure 
of the atmosphere. Likewise, the response of this line to 
the velocity gradient is also higher than in the other two. 
Then, a first step could be to characterize the response of 
all the lines to combined variations of the temperature, 
velocity and magnetic field and find observables (i.e., am- 
plitude ratios) that arc as insensitive to the temperature 
and velocity as possible and as sensitive to the magnetic 
field as possible. We are currently carrying out this study 
on realistic velocity fields, including shocks and temporal 
variations. 

The polarization amplification mechanism that we 
have discussed in this paper is not limited to plane- 
parallel atmospheres, although its effect is surely more 
important in plane-parallel atmospheres than in three- 
dimensional ones. The reason is that i) gradients in a 
three-dimensional atmosphere are weaker given the in- 
creased degrees of freedom and ii) non-resolved motions 
or large variations in velocity direction along the medium 
mix the contribution of different layers and broaden the 
profiles, diminishing them in amplitude. In any case, 
strong three-dimensional velocity gradients might cre- 
ate preferred directions along which the plasma becomes 
more optically thin through the radiative uncoupling 
mechanism discussed in this paper. 

The amplification of the radiation field anisotropy 
through vertical velocity gradients is a general and in- 
teresting phenomenon that improves our understanding 
of the stellar atmospheres. With the present investiga- 
tion we have obtained some feeling about the formation 
of the Ca II IR triplet in dynamic situations. 
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APPENDIX 
A. SPECIAL FUNCTIONS IN SECTION 3 
Introducing Equation \W\ into Equation [2] 
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Introducing the variables x = (y 1 — vq)/w, a — Avd/w, and £ = v z vq/(cw), then the mean intensity in the comoving 
frame 
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In passing from Equation (|A1[) to Equation (|A2[) . we have extended the integration limit on a; to oo. Analogously for 
the anisotropy in the comoving frame 
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where we have made use of 
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From them, the values for Jq and the anisotropy ^/2Jq/ J$ are trivially derived. 
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In the high velocity limit (£ —} oo), 2o(a;£) = (2 — u )/2, and Z 2 (ct;£) = u/4 (regardless of a). In the low velocity 
limit: 
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B. STATISTICAL EQUILIBRIUM EQUATIONS. 

The rate equations for the considered problem are as follow. They have been obtained by particularizing to the 
model-atom of Fig. 1 the equations contained in Sects. 7.2 and 7.13 of Landi Degl'Innocenti & Landolfi (2004). 
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where A u £ and Bi u are the Einstein emission and absortion coefficients; Ci u and C u e. ar c the excitation and dcexcitation 

(2) (2) 

inelastic collisional rates, respectively; Cg and C u ^ are the collisional transfer rates for alignment between polarizable 
levels (with J > 1/2); and D\ is the depolarization rate of the K-th multipole of level i due to elastic collisions with 
neutral hydrogen atoms. The p$ elements are referred to a coordinate system with the quantization axis along the 
solar local vertical direction. 

C. TWO-LEVEL ATOM CALCULATION IN A MOVING ATMOSPHERE. 

Figure CI is similar to Fig. [5J but for velocity fields with £ = 5 and maximum gradients occuring at different 
positions along the atmosphere. 
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Fig. 10. — Calculations in isothermal two-level atom moving atmospheres with £max = 5 . We assume a very strong line (r c = 0) and 
e = 10~ 4 . The highest velocity gradient occur at r; = to = 100, 10, 1,0.1,0.01 for a, b, c, d and e, respectively. The case f corresponds to 
the solution in a static atmosphere. The vertical dotted line marks the position of T\ = 1. 
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